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ABSTRACT 


The develc^nnent of pair-photon cascades initiated by high energy 
electrons above a pulsar polar cap is simulated numerically. The 
calculation uses the energy of the primary electron, the magnetic field 
strength, and the period of rotation as parameters and follows the 

curvature radiation emitted by the primary, the conversion of this 
radiation to e+ e- pairs in the intense fields and the quantized 
synchrotron radiation by the secondary pairs. A recursive technique allows 
the tracing of an indefinite number of generations using a Monte Carlo 
method. Ganina ray and pair spectra are calculated for cascades in 

different parts of the polar cap and with different acceleration models. 
We find that synchrotron radiation from secondary pairs makes an in^rtant 
contribution to the ganria ray spectrum above 25 MeV, and that the final 
gamma ray and pair spectra are insensitive to the height of the 
accelerating region, as long as the acceleration of the primary electrons 
is not limited by radiation reaction. 
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I. INTRDDUCTICN 


Pulsar magnetospheres are believed to be extremely favorable 

iz. 

environments for particle-photon cascades. Strang magnetic {>^ 10 gauss) 
and electric (<^10 esu) fields near neutron stars can provide all of the 
necessary ingredients for these cascades: high eneigy particles, efficient 
radiation processes, and a pair production mechanism. In the present 
picture, particles are accelerated by rotation- nduced electric fields 
above the polar cap and move along curved magnetic field lines. If the 
resulting curvature raxiiatioii from thes»> particles extends to gamma ray 
energies " 1 GeV, then the curvature photons have a high probability of 
producing electron-positron pairs in the magnetic and electric fields. The 
pairs resulting from magnetic pair production will in general be created 
with momentum components perpendicular to the field. Each member of the 
pair will therefore radiate synchrotron photons, many of which will be able 
to produce a second generation of pairs. This process will continue until 
the synchrotron photons are no longer energetic enough to pair produce. 
The surviving photons then constitute the observable high energy radiation 
from the pulsar, and the pairs may be input to a coherent process wnich 
produces the radio (and possibly optical) emission. 

The importance of such cascades in accounting for the observed radiation 
from pulsars was originally pointed out by Sturrock (1971) . Most current 
pulsar theories rely on electromagnetic showers for the production of the 
low energy particles which radiate coherently at radio wavelengths 
(Ruderman and Sutherland 1975, Cheng and Ruderman 1977, Arons <and 
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Scharlemann 1979) . Although these models differ as to how the accelerating 
electric fields are induced above the poles of the star, they all agree 
that pair production cascades will limit the potential drop below its 
maximum possible value. 

Most of the models which have been pressed to explain gamma ray 
emission from pulsars focus on the curvature radiation from the accelerated 
primary particles and the attenuation in the magnetosphere due to pair 
production. Hardee (1977) showed that particles accelerated by the polar 
gaps proposed by Ruderman and Sutherland (1975) (hereafter referred to as 
RS) were capable of producing curvature gamma rays with energies in excess 
of 1 GeV. A detailed numerical model involving primary curvature radiation 
and pair production in the magnetosphere has since produced calculated 
gamma ray spectra which give good fits to the observed spectra of the Crab 
and Vela pulsars (Harding, Tademaru and Esposito 1978; Harding 1981). 
These calculations ignored the details of the acceleration mechanism and 
the radiation produced while the primaries were being accelerated. But 
analytical calculations of the spectrum produced by particles in the 
accelerating region (Ayasli and Ogelman 1980, Ayasli 1981) have shown that 
radiation reaction during the acceleration cannot be important in fast 
pulsars like the Crab and Vela, or the spectrum would be much flatter than 
observed. Therefore, the gamma radiation from the acceleration region 
should not be an important contribution to the total spectrum. It seems 
that, as a first approximation, the observed gaiuna ray spectra of the Crab 
and Vela pulsars are best explained as curvature radiation from high energy 
primaries which are losing energy (ie. are above the acceleration region) . 
A calculation by Salvati and Massaro (1978), which made an attempt to 
include synchrotron radiation from the first generation of pairs suggested 
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that synchrotron gamma rays may also make an in^rtant contribution to the 
spectrum at energies below 1 GeV. 

In this paper, we present a scheme for numerically simulating the 
cascade produced by a primary electron, which we have used to calculate 
the resulting gamma ray and pair spectra. The technique allows a 
comparatively accurate treatment of geometrical aspects of the cascade 
process, to which the resulting spectra are extremely sensitive. The 
calculation follows the curvature radiation emitted by the primary, the 
conversion of this radiation to electron-positron pairs, and the subsequent 
quantized synchrotron radiation by the secondary pairs. A recursive 
procedure is aole to trace an ndefinite number of higher generations of 
pairs produced by the synchrotron ganma rays. We investigate the 
contribution of radiation from all generations of pairs to the gamma ray 
spectrum above 25 MeV, the number of pairs produced per pr inary electron, 
and the energy distribution of the pairs. Details of caf.cade development 
in the acceleration region and the angular distribution c>f pairs and gamma 
rays from the cascade will be treated in another paper. 


II. Numerical Simulation of Pulsar Cascades 

In the model considered here, the magnetic field is assumed to be purely 
dipolar above the stellar surface. Electrons are drawn from stellar 
surface layers at the polar cap zones, and accelerated along the magnetic 
field lines up to a height ccxjparable to the polar cap radius. Above this 
height, the magnetospher ic plasma is assumed to cut off further 
acceleration, and the electric fields induced by the co-rotating plasma 
have no significant parallel component. 



” 6 “ 


ihe calculation described here traces the cascade initiated by a single 
primary electron which or ginates from a given point on the polar cap. Its 
input parameters include the pulsar period and radius ^ the surface field 
strer^th and orientation of the magnetic field relative to the rotation£d 
axis, and the peak .energy achieved throiigh acceleration. IVk> versions of 
the calculation have been developed in order to treat both intantaneous and 
continuous acceleration. 

A full cascade resulting from large numbers of primary electrons, whose 
energetic and spatial distributions over the polar cap zones are supplied 
as additional input parameters, can in principle be found by taking 
weighted results from separate single particle calculations of the type 
discussed here. Unfortunately, while the solutions found by Deutsch 
(1955) for the surface field configurations of magnetic stars in vacuo 
might suggest approximate forms for these distributions, ti*e surface charge 
structure in pulsars surrounded by magnetospheres is not well known. 
The spectral results reported in Section III below are thus restricted to 
the case of primary acceleration from the outer rim of the polar cap. 
However, it will be shown that this case may be a good approximation to the 
actual time integrated spectrum from the polar cap. 

Ihe calculation involves two different reference frames; the rotating 
frame where the star and dipole field are stationary, and the nonrotating 
or observer's frame. Tracing of particle trajectories takes place in the 
rotating frame because they follow field lines, while tracing of photons 
takes place in the nonrotating frame where they follow straight line 
trajectories. 
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a) Primary electron trajectory 

In the first version of the calculation, the enerqy of the primary 
electron at a given height above the stellar surtaK:e is tipecified eus input. 
Above this height, the electron is assumed to "coast” outward along the 
field line corresponding to the chosen surface coordinates, losing energy 
via curvature radiation. The second version, which considers a continuous 
acceleration of the pri.nary electron to a height comparable to the polar 
gap heights suggested by PS, has also been studied. It turns out that as 
long as the acceleration region does not extend much further above the 
polar caps, the final ganma cay spectra which result from primary electrons 
with similar peak energies is not very different in the two versions (see 
Section III) . 

In both versions of the calculation, the outward trajectory of the 
primary electron is divided into segments of arc length as, over each of 
which the electron loses (at most) a fixed fraction ■ dE / E of its 
present energy, E ■ 7 me”, due to curvature radiation. If the path length 
over which this fractioned energy loss occurs should become large compared 
to dimensions over which the fields change, the segment length is 
limited to a fixed maximum value as^^^, measured in units of the stellar 
radius R,. In the version which includes acceleration of the primary 
electron, s^^^ must be further limited to an a^Ptepriate fraction of the 
acceleration height. The value »| > 0.1 has been found sufficiently snail 
to be adopted as a standard vedue in the calculations reported here. 
Values adopted for As are listed in Tedole 1(a). 
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The division into path segments is made by finding the exact radius of 
curvature , 


( 




K (. ' *■ 





( I V 
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for the dipole field line, r ■ K sin C- , at the electron position (the 
starting point of the segment), and then estimating the total energy loss 
rate. 



due to curvature radiation at that point (see for exarple Jackson 1975), 
where r, * e /m c is the classical electron radius. From this rate, the 
increment of arc length over vAiich (d Y / dt) (ds /c) “ 

derived. (If the calculated increment d s is larger than the maximum value 
this maximum becomes the increment and the specified energy loss is 
replaced by the smadler loss over this increment.) 


Ihe energy loss from curvature radiation is subtracted and the gain due 
to acceleration, if any, is added to the energy at the start of the 
segment, so that the net energy change is 
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E is the electric field, parallel to the magnetic field, which 

’ I 

accelerates the electron. We assume the simplest model where is 
constant, 

I 1 u 
> u 


where V is the potential dr<^ and h is tt» extent of the acceleration 
region above the stellar surface. In the first version we consider, V ■ 0 
(so that 0 everywhere) and the electron trajectory starts at the 
stellar surface with initial erwrgy E^. equsd to the maximum energy. In the 
second version, values are adopted for V and h, and the electron starts at 
the stellar surface with E. « m c^ (ie. at test), where it accelerates from 

V 

s » 0 to s ■ h along its msignetic field line. If the first term remains 
small cotppared to the second term in equation (3), then the acceleration 
will not be radiation reaction limited (NRR case) and the electron's 
maximum energy will be E^ ■ V. If, however, the first term becomes 
conparable to the second term, the acceleration is radiation reaction 
limited (RR case) , and the electron will achieve a itaximum steady state 
energy. 




less than V, in a distance. 
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After equation (3) is used to find the new energy, the position is shifted 
along the dipoie field line by the distance ^s. For this purpose, a table 
relating arc length displacenents from the stelJar surface to coordinates 
(r, h) along the chosen field line is established at the beginning of the 
calculation. 


b) Curvature radiation 

Before the calculation proceeds to the next segment, the outward 
pr(^>agation of the curvature radiation emitted over the present segment, 
and the ensuing shower of secondary photons and electrons which it 
generates, are completely treced. This task is accomplished by a recursive 
procedure, whose inputs are the energy, initial position and time 
coordinates, and the direction cosines of a photon emitted at any given 
stage in the cascade development. Figure 1 schematically illustrates the 
cascade produced by a curvature radiation photon emitted during the first 
step in a primary electron trajectory. 

At its lowest level, the procedure fellows the "burxile" of curvature 
pht'tons emitted from each segment. Since the segments over which most of 
the curvature radiation is emitted are quite short compared to 
characteristic lengths such as or the radius of the polar cap, the 
calculation assumes that all photons are emitted from the midpoint of the 
segment. Furthermore, the direction of emission (in the rotating frame of 
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the magnetosphere) for the entire ^ectrum of photons is taken to be that 
of the tangent to the field line at this point (since '^^1/-^ « 1). in 
the non-rotating frame, \^ere the calculation follows the photon to its 
pair conversion point, there is an aberrati(xi of the enissicxi direction so 
that the photon will not initially be propagating tangent to the magnetic 
field in that frame. 


The spectrum of curvature radiation may be found from classical 
electrodynamics (it is in fact identical to that of classical synchrotron 
radiation froir an electron with the same radius of curvafire) and may be 
written in the form (Jackson 1975) , 


c\T.. 
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is the cutoff frequency, and the function K(x) is defined as, 
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It has the asin^ptotic forms: 
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which we use to evaluate the function in these limits. Around x » 1/ the 
function values are interpolated from a table calculated values. This 
spectrum is divided into increments A*-j, and the "number" of photons 
in each spectral increment is estimated by dividing the energy, 
(Jt/. I ' contained within the incr«nent by the average photon er«rgy 
of that increment. The A photons are all assumed to have this average 
energy. Standard values for the ^xsctral parameters chosen in this 
calcuJ.ation are given in Table 1(b). 

The calculation in effect takes only one '“f the photons in each 

spectral increment and traces the cascade which that photon produces. Pair 
conversion is assumed to occur when its probability of survival,? (x) * 

exp(-/(dx /A )], after propagating a distance x through the anhient 
fields, has dropped to 1/e. Here, A is the local mean free path for pair 
production. However, it is assumed that the cascades resulting from the 
other photons emitted within this same spectral inc: snent at the same point 
are identical. lhat is, the secondary photon and electron spectra 
resulting from the one sample photon are multiplied by An^^ before adding 
their contribution to the accumolated gamma ray spectrum. 

This use of the relative numbers of photons in each increment as a 
weighting factor is done only for the curvature radiation photons, since 
the synchrotron radiation emitted subsequently by secondary electrons is 
both qiiartized and random in azimuthal direction, and requires the tracing 


| A C ^ = \ 

I 
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of individual photcns. The assunption that pair production occurs vhen 
whai P(x)=l/e is therefore made only for the curvature photons. Padr 
production points for the synchrotron photons are determined by giving 
them each a random probability of survival. A pure Monte Carlo approadi 
vioiild also require the tracing of curvature photons on an individual basis , 
but this averaging approxiination has been adopted to a^roid prohibitive time 
requirements. The principal justification for this approach lies in the 
fact that the mean free paths for the energetic curvature photons are 
typically short coopared to the stellar dimensions ^ so that the dispersion 
in these distances may be neglected in tracing the cascade development. In 
addition r it may be noted that the final cascades produced bv curvature 
photons are typically quite large, and hence those cascades initiated by 
parent photons which are identical in energy, position, and direction 
should be statistically similar. 

c)- Pair production 

We have investigated several iiethods for following the outward 

propagation of the photons through the magnetosphere. This procedure is 
quite sensitive to the choice of path-length increments, since the 
attenuation coefficient for magnetic pair production is a strongly varying 
function of the strength and orientation of both the magnetic and electric 
fields. In addition, the photon follows essentially a straight-line 

trajectory in the nonrotating frame (ignoring general relativistic effects 
and slight deviations due to refractive properties of the intense field) , 
wh''e the star and the near -zone fields rotate beneath it. As mentioned 
earlier, the aberra*-ion of the emission direction must also be taken into 
account. The inclusion of photon aberration is crucial in correctly 
determining the initial pair production rate (Harding et. al. 1978) . All 
of these effects must be treated carefully in order to arrive at reasonably 
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accurate pair conversion points. Of the methods explored to date, the most 
satisfactory approach has been to use initially small increments in the 
photon path lengths , which are dynamically rescaled to greater (and less 
time consuming) values as permitted by the step-by-step changes in the 
ambient fields. Ihe parameters adopted in this procedure are shown in 
Table 1(c). 


The expressions used here for the pair production attenuation 

coefficients are those found by Daugherty and Lerche (1975), which include 

the modifications to the attenuation coefficients in pure magnetic fields 

(Toll 1952; Erber 1966; Tsai and ErL«r 1974) due to the presence of 

z - 

perpend ic'jlar electric fiel^Js fwtere B - E >0 and E-B ■ 0). Without 
losF of generality, a frame may be chosen such that B » B z and E ■ E y, in 
which case the (unpolarized) attenuation coefficient for a photon with 
arbitrary direction cosines ( »'j ^ ^ ^ ^ /| ) has the asynptotic form. 






I V 







(■7 ) 


in the limits B/B << 1, E /me >> 1, where B * m"c*/e ft is the critical 
field strength. 


Although the surface acceleration model used here implies that the above 
condition E-B « 0 is not satisfied everywhere, the validity of equation (7) 
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is not significantly affected even when E^B does not vanish. In fact, the 
effects of parallel components of eleccric fields have been found to be 
negligible for the photon optical depths (Harding et. al. 1978; Ayasli 
1978) , evai f'?r the upper limits based on the vacuum model (Deutsch 1955) . 

At the conversion point, both mendsers of the pair are ass.gned the same 
direction and half the energy of the parent photon. These high energy 
approximations are reasonable (Ayasli 1978) for the cascade model, but have 
been made here for the equally cogent reason that no detailed investigation 
of the angulcir and spatial distributions in magnetic pair production is 
known to us. If the present calculation should be extended to moderate or 
low energy gamma ray emission, such information would probably be required 
to derive accurate spectra. 

d) Quantized synchrotron radiation 

The initial pair direction is used to estimate the initial particle 
energy tranverse to the local magnetic field. In general, this energy will 
be considerably higher than the local (quantized) ground state energy, and 
each member of the pair must rapidly lose this energy through a burst of 
quantized synchrotron radiation. It is iin»rtant to note that, in ccxitrast 
to curvature radiation, the photon energies from synchrotron radiation in 
superstrong magnetic fields are typically coitparable to the energies of the 
particles themselves. Moreover, the emission rates are sc rapid (even 
after the Lorentz time-dilation effects due to parallel motion are taken 
into account) that essentially all the initial tranverse energies of the 
created pair are lost before the particles can move a significant distance 
outward along the local field line. For both of these reasons, the 
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contribution of the synchrotron c<»tponent to the emergent gamma ray 
spectrum has been found to be quite significant. (In contrast, the 
curvature radiation from all generations of the seconda*./ electrons has 
been estimated to make a negligible contribution to the final spectra, at 
least above 25 MeV. p\or B ^ 8 X 10^^ gauss, syncdirotixMi radiation from the 
pairs ^dway8 dominates, because V'sin4' y 1 for the photons which pair 
produce.) An additional consequence of the rapid emission rates is that 
acceleration of the pairs during their synchrotron burst will be negligible. 

In order to follow each manber of the pair through its synchrotron 
emission sequence, the equation of the dipole field line which passes 
through the ccxiversion point is found and the particle momentum parallel to 
the field line is inferred from its initial angle to the field, V . This 
parallel momentum in fact determines a "longitudinal" Lotentz factor 

y - { 4 — - CCS ^ ^ 

<1 ^ y “ 

i 

which is a conserved quantity under the aj^toxiitiatiors adopted here. For 

^'This calculation hats assumed that the synchrotron photons are aill 
emitted in the orbital plane of the electron, in the frame where the 
electron moves transverse to the B field. In this ^prcxijnation , is 
strictly conserved since the recoil of the electron when it emits a photon 
cannot change its energy paurallel to the field. We have verified numerically 
that is conserved in this case. 

M 

each raerber of the pair in turn, the caLLculation then follows the quantized 
syrichrotron emission sequence as the electron drops toward its ground state 
orbital in the local field. For each photon emission, a series of Lorentz 
transformations are made to orient the coordinate axes along the field line 
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aiul remove the longitudinal motion, emit a photon in the frame of purely 
"circular" motion, and then return the energy and direction of the emitted 
photon in the original (stellar) frame. I^e quantized emission spectrum 
given by Erber (1966), 



r 


has been used to sanple the photon energies in the transverse frame. The 
K (x) function [cf . eqn. (6) ] is the same as that ag^aring in the 
curvature radiation spectrum. The photons are assumed to be emitted in the 
plane of the electron orbit with random azimuthal angles. 


The recursive aspect of the calculation appears at this point, since for 
each anitted synchrotron photon the entire [procedure of tracing the photon 
propagation to its annihilation point (or until it has been determined to 
have escaped conversion), followed by the tracing of its resultant pairs 
through the synchrotron emission sequence, is begun anew. The recursion 
may be pushed to an indefinite number of levels, depending on the number of 
successive photon and pair generations emerging in the cascade. The 
process is eventually brought to a conclusion when the Ic^r energy photons 
belonging to the succeeding generations escape from the magnetosphere 
without pair producing. In addition to the energy degradaticn, the rapid 
decrease in field strengths with distance from the s:ar is also 
instrumental in ending the cascade. Table 2 indicates typical levels of 
recursion (generation numbers) encountered in calculations with varying 



surface field strengths and primary electron energies. 


After the calculation finishes tracing the conplete shower due to the 
emission of a given synchrotron photon# the transverse energy of the 
electron which produced that photon is reduced by the amount of the emitted 
photon energy# and the electron is displaced outward along its guiding 
field line by a distance estimated from the toted synchrotron emissicm 
rate# 


A X 




Here# G(r) is the loss rate function defined by Erber(1966). 

The synchrotron emission process is then repeated for this electron 
until its energy drops below a suitable threshold# taken here to be 100 
MeV. When this process has been completed for each member of the pair in 
turn, the procedure returns to the next lower level of recursion. At the 
lowest level# namely that of the primary electron, the electron is 
displaced along its outward trajectory as described above, and the next 
bundle of curvature photons is traced. 

The calculation accumulates a spectral distribution of all ganma rays 
above a certain cutoff energy (here taken to be 25 MeV) which are found to 
escape pair conversion in the magnetosphere. Spectral distributions are 
also obtained for the accumulated ger«rations of pairs, in which the 
"final" energy of each electron or positron is taken to be its energy 
parallel to the magnetic field, , which is conserved. Ttie final energy 
of pairs produced in the accelerating region is# 
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SO that electrons (+sign) are accelerated way from the star and positrons 
(-sign) are accelerated back towards the stellar surface. This gain or 
loss of pair energy parallel to the field will not affect the sydirotron 

radiation, because of the short radiation timescales, but if a member of the 
pair should gain enough energy, its curvature radiation may contribute to the 
fined spectrum. We do not consider this contribution in the presait version 
of the ccdculation. Ihese spectra are found mainly as a check on the overall 
energy distributions among particles in the final cascade. No attempt has be^ 

made here to trace their subsequent (and much lower) losses due to 
curvature radiation, nor to guess at any subsequent acceleration these 
particles may undergo as they approach the light cylinder. In addition to 
the spectral distributions, the multiplicity, M «{N^ + )/N ^ , or the 

ratio of secondary to primary particles, is also found. 

III. RESULTS 
a) Gaima-ray spectra 

Ihe first version of the calculation, as described in Section II, 
ignores the details of che acceleration and assumes that the primary 
electron is instantaneously accelerated to energy E^at the stellar surface. 

As we will show later in this section, acceleration over a finite distance 
for parameters of current polar cap models does not change the calculated 
spectrum significantly. Figure 2 shows the build-up of the primary 
curvature spectrum and the accumulation of escaping photons from the 
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cascade as the primary electron steps out along its trajectory. The 
curvature spectrum for the first few steps has the characteristic -2/3 
slope but, as the primary loses energy, the shift of the curvature cutoff 
to lower energies steepens the spectrum to an asyn^Jtotic slope of -5/3. 
Because the radius of curvature of the field line also increases as the 
electron moves outward, the contribution to the final curvature reidiation 
spectrum from successive steps decreases. We find that for the parameters 
listed in Table 1, the contributions beccmie negligibly snail above the 20th 
step. The cascade gamna ray ^ctra from the first 3 steps have extremely 
sharp high energy cutoffs, indicating that virtually all photons above 500 
MeV (and many below 500 MeV) are efficiently converted to pairs and lower 
energy photons. In fact, in these first few steps near the stellar 
surface, the gamma ray spactrum is almost entirely dominated by synchrotron 
photons. Curvature radiation from higher steps farther away from the star, 
where the fields have decreased and photons can escape more easily, fills 
in the high energy part of the final spectrum. The final high energy 
cutoff is a balance between the decreasing optiCcLL depth of the 
magneto^hete (the increasing energy of the escaping photons) and the 
decreasing cutoff energy of the curvature radiation along the trajectory. 

The calculation shown in Figure 2 followed the cascade initiated by an 
electron moving on the outermost open field line (ie. at the edge of the 
polar cap) . Figure 3 shows the final gamma ray spectrum from 10 eV 
electrons on field lines originating at different magnetic colatitudes, , 
on the stellar surface [in units of the polar cap half angle defined by 
sin » (R .-/c)^^]. There are basically two effects whicn change the 
spectrum as a function of , and they are both due to the magnetic field 
line radius of curvature. First, the increase in radius of curvature 
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toward the magnetic pole gives a lower curvature radiation cuto££i so there 
ctte fewer high energy photons to feed the cascade. Second, the larger 
radius of curvature lowers the pair production optical depths, since 
photons will have to travel further to encounter a sizable perpendicular 
field coniionent. This allows higher energy photons to escape, decreasing 
the synchrotron contribution to the final spectrum. Clearly, the most 
extensive cascades are produced by electrons accelerated near the edge of 
the polar cap. The observed time integrated spectrum from the pulseir would 
be an average of the spectra produced by electrons at different 6^ across 
the pulse, weighted by the primary flux across the polar ceqp. The case 

9 ■ is probably the best aE^roxiroation to the time-integrated 
spectrum since i) most of the gamma ray flux is produced in cascades from 
primaries with (>[, ■ ii) most of the area swept out by random observer 
lines of sight to a pulsar will be in the outer parts of the polar cap 
iii) in current models, acceleration is strcxigest at the edges of the polar 
cap iv) the variation of the spectral index with 9^ is small except near 
the poles. 

The dependence of the c^dculated gamma ray spectrum on pulsar period is 
shown in Figure 4. The steepening of the spectrum with decreasing period 
results primarily from the increasing magnetic colatitude of the outermost 
open field line (cf. Fig. 3) Also iirportant for short periods is the 
increase in the pair production optical depths due to the rotation-induced 
electric field. Because of these effects and the arguments given in the 
preceding paragraph, short period pulsars should have steeper ganvna ray 


spectra. 
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In the second version of the calculation, the primary electrons are 
accelerated to a ^^ecified energy over a finite distance above the polar 
cap. The contribution to the gamna ray spectrum from this accelerating 
region is shown in Figure 5 for different values of the voltage and the 
acceleration distance. Cases A1 and A2 have parameter vedues which are 
approximately those estimated by RS and Arons and Scharlemann (1979) 
respectively for the Crab pulsar. In each of these cases the accelerating 
particles do not experience strong radiation reaction ( the second term in 
Eqn. [3] dominates) and they gain an energy equal to the potential drop 
V. Case A3 is radiation reaction limited ( the first term becomes equal to 
the second in Eqn. (3] ) and the electrcwi reaches its steady state energy, 

t o 

about 10 eV after traversing only half of the accelerating region. The 
gamma ray flux from acceleration regions Al and A2 is negligible at all 
energies compeured to the gamna ray flux produced above the acceleration 
region. This is because, for values of h less than or conparable to the 
polar cap radius, the distance traveled during acceleration is much smaller 
than the typical length scale over which the particle loses its energy from 
curvature radiation. Furthermore, when the acceleration is not radiation 
reaction limited, the particle travels a relatively short distance near its 
maximum energy while it is being accelerated. For larger acceleration 
heights as in case A3, the electron spends more time at or near its maximum 
(steady state) energy and therefore, the contribution to the gamma ray 
specuum is greater. However, acceleration heights this large at the polar 
cap will be unstable against pair production breakdown ( RS ) . In fact, in 
case A3 we find pairs produced well within the region s < h. Polar 
acceleration regions will generally be non-radiation reaction limited, 
except possibly in long period pulsars. The pairs and gamma rays produced 
in the acclerating region are therefore not an important contribution to 
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the final g^uma ray or pair spectra. Hovwver^ they are undoubtedly 
iirportant in determining the structure and dynamics of the accelerating 
region, and we hope to consider this problem in future work. 

Figures 6a and ^ show gamma ray spectra for periods of the Crab and 
Vela pulsars, each calculated for three different values of the surface 
magnetic field strength. In these calculations, the electron is started 

t % 

from the stellar surface with 10 eV and does not experience any 
further acceleration alcx^ its trajectory. The effect of increasing the 
magnetic field is to steepen the spectrum by producing a more extemsive 
cascade. A higher field increases the pair production optical depths, 
ranoving high energy photons from the primary spectrum and converting them, 
through syrichrotron radiation, to low a\ergy photons. A shorter period 
amplifies this effect, as discussed previously, through higher pair 
production and decreased radii of curvature of the field. The observed 

I w 

COS-B data points for the Crab and Vela seem to best fit the 10 gauss 
spectra for both pulsars. 


b) Pair ^>ectra 

The calculation also keeps track of the total number of pairs producjed 
in a cascade as well as the number of pairs in each energy bin. Figure 7 
shows the differential number spectrum of secondary electrons and positrons 
from a cascade having parameters of the Crab pulsar. These energies are 
the computed energy components of each member of the pair parallel to the 
magnetic field at the tijne they are produced. Changes in the spectrum due 
to Compton scattering or annihilation may occur, but the present 
calculation does not include such effects. Subsequent radiation losses. 
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howevec, should be negligible unless the particles are accelerated. Ihe 
spectrum is roughly a power law with a cutoff at around 5 GeV. This high 
energy cutoff moves up for stronger fields, because the highest energy 
photons produce pairs at snaller angles , giving the pairs larger parallel 
energies. The low energy turnover is a real effect and not a result of the 
photon energy cutoffs. It results rather from the geometry of the fields, 
which determine very critically the lowest energies of the pair-producing 
photons and the angle to the field at which the pair is produced. The pair 
spectra for cascades in stronger fields have lower turno';er energies 
because lower energy photons are able to produce pairs. 

Results for the nniltiplicity, or the total number of electrons and 
positrons produced per primary electron in the cascade, for values of 
magnetic field strength ano period are shown in Figure 8. No pairs are 

IC 

produced for fields below 10 gauss for the period of the Crab. Above this 
field strength, the nunber of pairs rises sharply and then begins 

I 2, 

levelling off around a field strengtii of 10 gauss. The cascade 

multiplicities for the Vela period, P • .089 s, are lower by a factor 

of 2 or 3. Multiplicities edso decrease as the magnetic colatitude of the 

primary electron decreases for reasons discussed in Section Ilia. For 

I- 3 

typical field strengths around 10 gauss, multiplicities of 10 are easily 
achieved. While these numbers are twjch smaller than those originedly 
estimated for pulsar cascades by Sturtock (1971), they are in good 
agreement with those estimated by RS. The total energy Ln photons emerging 

I . 

frcxn the cascade exceeds the total energy in pairs for fields below 5 X 10 
gauss, but above this field strength the energy in pairs exceeds the energy 
in photons. This suggests that short period pulsars with high field 

strengths may lose a greater fraction of their rotational energy in 
particles than in radiation. 
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It is possible to estimate the primary and secondary particle fluxes 

from the Crab and Vela pulsars from these results r using the observed ganna 

u. 

ray fluxes to normalize the conpated gamma ray spectra at 10 gauss. 

Ihe primary electron flux required to give the observed gamma ray flux 
-( 

levels is 6 X 10 j for the Crab (assuming a distance of 2 kpc.) and 2 X 

? .| 

10 8 for Vela (assuming a distance of 0.5 kpc.). Iheae fluxes are of the 
same order of magnitude as the flux of particles expected from breakdown of 
a polar gap, which is 


_ 2 _ 




TT V' 




- I ^ 


0 


where N is the number of secondaries produced per primetry in the ge^ which 

are accelerated to energies comparable to the primary, and ■ 7 X 
IC -I -3 

10 P . an is the corotation charge density at the stellar surface. With 

the calculated multiplicities for a complete cascade, the secondary 

77 .| 

particle flux from the Crab would be 3 X 10 s and the secondary particle 

flux from Vela would be 4 X 10 s . This particle flux from the Crab is 

4 / -I 

almost 4 orders of magnitude less than the 10 s estimated by Sturrock 
(1971). 


IV. DISCUSSICN 

This paper has developed a technique for tracing electromagnetic 
cascades near pulsar polar caps. We have presented results of numerical 
calculations which predict tiie ganma ray and pair cpectra for a standard 
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polar cap model and have shown that the predicted ganina ray spectra for tlie 
Crab and Vela pulsars are in good agreement with the observed spectra. 
Furthermore / it has been shown that a polar gap origin for the observed 
ganana rays less tivui 5 GeV is energetically feasible. This conclusion 
differs from that of Ayasli and Ogelman (1980) 'tio ocmclude that polar 
gaps or acceleration regions near the poles make a negl^oible contribution 
to the total ganr* ray emission conparec to the outer gaps which have been 
proposed by Cheng and Ruderman (1977) . They argue that the number of pairs 
involved in eadi polar gap breakdown is a factor of 10 P lower tium tl'ie 
number of pairs from an outer gap breakdown. Since these pairs would 
produce much of the ganna ray emission around 100 NeV, the ganma ray 
mission from the polar gap would, in thoir model, be lower tiian that from 
the outer gap by the same factor. However, they neglect the contribution 
from secondary pairs produced by cascades from the accelerated prinv*ries 
above the polar gap. From our results, this synchrotron contribution 
provades about a factor of 10 increase in conversion efficiency of primary 
electron energy to ganina rays in the Crab, which coR|:er.^ates for the lower 
particle flux from the poleir gap. 

W have found that the ganna rays produced while the primary electron is 
accelerating to its maximum energy make a negligiWe contributiem to the 
final gazma ray spectnxn for stable polar accelerat:ing potentials. Higher 
potential drops or lerger acceleration regions would make more svibstantial 
cont::;ibutionc, but would probably be unstable. Exactly v*iere and hew these 
accelerating potentieds break down are problans which could be 
investigated in detail with the technique developed here. A particularly 
interesting question is whether the potential drops are limited by 
synchrotron - pair production cascades at the top of the e<cceleration 
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cegion, thus producing a quasi-steady discharge as suggested by Arons and 
Scharlemann (1979 ) , or whether the curvature radiation - pair production 
cascades proposed by RS can completely short out the parallel electric 
field. In the first case, the voltage and particle fluxes are constant in 
time while in the second case, they are time variable. Which case actually 
occurs in pulsars is of central iirportance to radio mission theories, and 
it hinges on whether the secondary pairs can be accelerated to primary 
energies and radiate high energy curvature photons. 

The presence of higher magnetic field nultipoles has been suggested in 
pulsar nodels. These nultipoles could affect the near field structure by 
decreasing the radius of curvature of field lines. This would cause an 
increase in both the energy of the curvature photcxis and in the pair 
productiOTi rates in cascades near the surface, possibly producing steeper 
spectra than we have v^culated for pure dipoles, and certainly producing 
nore extensive cascades. 

The result that cascades at different magnetic colatitudes produce 
different gamma ray spectra suggests that there should be observable 
spectral variations across the ganma ray pulse. There is observational 
evidence for this in the Vela pulsar (Kanbach et. al. 1981), although it is 
not clear where the gaima ray pulses are originating because of their phase 
lags relative to the radio pulse. Detailed predictions of spectral 
variations with pulse phase would require a more complete model of cascade 
development across the polar cap. The variations found in the flux of 
gamma rays with magnetic colatitude also suggests that the possible 
observed flux variations in the Crab pulsar (Greison et. al, 1975; Wills 
1980) could be due to neutron star precession, presenting different parts 
of the polar cap to the observer at different times. 
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The calculation described in this paper could be refined and extended, 
in particular to cover lower gainna ray energies and to investigate pulse 
shapes and ternporal behavior. Those approximations which have been made to 
avoid long con|Xitation times may be refined without great difficulties. 
However, several of the approximations mentioned above have been necessary 
because the physical processes involved are not sufficiently well 
understood. In the case of synchrotron radiation in super strong magnetic 
fields, convenient analytic forms for the fully quantum mechanical regime 
are not available. The expressions used in the present calculation are 
inaccurate when the energy of the emitted photon approaches the energy of 

the electron, a situation quite probable for field strengths in excess of 

12 - 

10 gauss. The approximation that the energy of the electron parallel to 
the magnetic field is conserved may also be incorrect when radiation 
reaction effects are important (Shen 1979) . In this case, the angular 
distribution of the emitted photons must be considered, so that the recoil 
of the electron parallel to the field may be taken into account. In the 
case of pair production, the energy and angular distribution of the pairs 
is not well known. The high energy ^roxurations made in this calculation 
(that the pair share the energy of the photon equally and that they 
initially travel in the same direction as the parent photon) are good only 
for photon energies E , » m c " . These fundamental processes must be 
investigated in greater detail before the technique described here can be 
extended to lower ganma ray energies. 


We would like to thank D. J. Thompson and S. Ayasli for discussions and 


helpful comments. 
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TABLE 1 

CASCADE PARAMETERS 


a) Primary electron trajectory 

Number o£ steps 20 

Maximnn fractional energy loss per step 0.1 

Maximum st^ length (As. < h: RR s VS 

’ a NRR h/5 

[^>* h R-. 


b) Curvature spectrum 


High frequency cutoff ( i /uJ^.) 7.0 

Low frequency cutoff ( a ) 0 

Spectral increments per dec^e 10 

Absolute low energy cutoff 100 KeV 


c) Pair production 


Initial step length in photon path 1 cm 

Factor by which steps increase 10 

Low field cutoff (B / 10"^ 

c»- 


d) SynchrotrOTi spectrum 


High frequency cutoff E,- 

Low frequency cutoff 25 MeV 

Spectral increments per decade 3 


* 
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TABLE 2 

CASCADE PHOTON GENERATIONS^ 


INITIAL PRIMARY MAGNETIC FIELD STRENGTH 

ENERGY B (GAUSS) 

E^, (eV) 



5 X 10 



5 X 10 

1 X 10'* 

2 X 10* 
5 X 10‘* 



12 . 

5 X 10 

4 

4 

7 
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figure CAFTICNS 


Figure 1 - Schematic representation of a pair photon cascade 
initiated by an accelerated electron above the 
polar cap of a pulsar. Four generations of photons, 
including both curvature radiation photons, and 
synchrotron radiation photons, are shown. 

Figure 2 - Cumulative spectra of primary curvature radiation 
photons (dashed line) and cascade photons (solid 
line) for various steps in the primary electron 
trajectory. 

Figure 3 - Final gamma-ray spectra from cascades initiated 

by 10 cV electrons at different magnetic colatitudes, 
9p in units of the polar cap half-angle, for a 

I L 

surface field of 10 gauss and period of .033 s. 
Figure 4 - Final gamma-ray spectra for different pulsar periods 

i 

cdl with initial electron energy 10 eV, surface 

\Z- 

field 10 gauss and magnetic colatitude 6 * 1.0 

Figure 5 - Gaitna-ray spectra irom the accelerating phase 
of the primary electron trajectory compared 


f 
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to the final gansna-ray spectrum produced in 

the case of instantaneous acceleration. Case A1 

I 'h 

heus V * 10 eV and h ■ 10 cm., case A2 has 

( > Cf 

V « 10 eV and h « 10 an., case A3 has V ■ 

1 5 ‘j 

2 X 10 eV and h * 10 cm, while case A0 

i 3 

corresponds to V « 0, h > 0, and E.« 10 eV. 

Figure 6 - Ganma-ray ^»ctra for the periods of the 

a) Crab and b) Vela pulsars and three different 

values of the surface field strength in units 

tz. 

of 10 gauss. Also plotted are observed spectra 
for each pulsar from Lichti et. ad. (1980). 

Figure 7 - Differential e+ e- spectrum from a cascade with 

i5 

initial electron energy 10 eV, surface field 

1 1 ,, 

10 gauss, magnetic colatitude ‘J . « 1.0, and 

period of the Crab, P » .033 s. 

Figure 8 - Cascade multiplicity versus surface field strength 
for the periods of the Crab and Vela pulsars and 
also for a more typical pulsar period of 0.5 s. 

I i •* 

All calculations have £. = 10 eV, and ‘y = 1.0. 
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